Gluons at finite temperature in Landau gauge Yang Mills theory 
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1. Introduction 

One of the most intriguing properties and challenging problems in QCD is confine- 
ment. In r ecen t years, the understanding of gluon confinement in Landau gauge has 
progressecP^l, j n a combined effort of methods using Dyson-Schwinger equations 
(DSEs)^ renormalization group (RG) met ho d and in a multitude of lattice 
calculations (for the most recent results, see e.g. Bowman et al. and Oliveira et 
al®). The question of quark confinement, on the other hand, is not yet at the same 
staged However, significant evidence exist s tha t gluon confinement is not affected 
by a small (physical) number of light quarkJ^^. Furthermore the non-perturbative 
features of QCD are most probably generated in the gauge sector. It is therefore 
reasonable to use the knowledge on gluon confinement also in applications, like the 
investigation of Yang-Mills theory at finite temperature. 

Several decades ago it was arguedUl that QCD and Yang-Mills theory undergo 
a deconfining phase transition at sufficiently high temperatures. This prediction 
has been investigated in many experiments and in lattice calculations . One of 
the striking expectations was that the high-temperature phase should consist of 
weakly interacting quasi-particles. However, although thermodynamic bulk proper- 
ties seem to logarithmically approach this behavior I I, the microscopic properties, 
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especially in the chromomagnetic sector, do not seem to behave as expected-^. Also 
there is not yet a theoretical explanation which can consistently explain all exper- 
imental obs ervations . Many attempts to reconcile this problem, e.g. hard thermal 
loops 

(HTL jnniTB] 

have been made, but did not succeed. This was to some extent 
anticipated, as the infinite-temperature limit of Yang-Mills theory is eq uivalent to 
a 3-dimensional (3d) Yang-Mills theory coupled to a massive Higgs fielcPEni 3d 
Yang-Mills theory is confining^* 1 , and thus non-perturbative interactions have to 
be expected. 

Here, recent results on this subject are reviewed, especially from DSEsPHIEl 
With a direct handle on gluon confinement available, it is possible to trace the 
fate of confinement along the complete temp era ture axis. As expected gluon con- 
finement prevails below the phase transitiorCpl In the high-temperature phase, 
part of the gluon spectrum is still confined^. This confinement persists, as an- 
ticipated, in the i nfinite- temperature limil^l These results agree with findings in 
lattice calculation :) 25 ! ^ ! 

This review is organized as follows. In section |21 basic properties of Yang-Mills 
theory are collected, including signals of confinement. Aspects of the formulation 
at finite temperature are presented in section [21 together with a discussion of the 
DSEs. In section 0] the results are given. These are discussed in section [SJ which 
includes the impact on the phase diagram and the possible influence of quarks. 
This section also concludes the review. 



2. Aspects of strong Yang Mills theory 

The theory studied is an equilibrium Yang Mills theory, governed in the Matsubara 

11 ly*7l 

formalism by the Euclidean Lagrangiar i 1 

£ = + c a d»Dfc b (1) 

F%, = dp A" - duA% - gf ahc AlA c v 
Df = S ab d^ + gf abc Al . 

Here F" denotes the field strength tensor, D ab the covariant derivative, g the 
gauge coupling, and f abc the structure constants of the gauge group. A® is the 
gluon field and c a and c a are the Faddeev-Popov ghost fields, describing part of 
the intermediate states of the gluon field a . The gauge chosen is the Landau gauge 
for technical reaso ns, most importantly because of the non-renormalization of the 
ghost-gluon vertejl^l At finite temperature no comparable set of investigations 
in other gauges have yet been completed. There are i nves tigations in the vacuum 
in Coulomb gauge using DSEs 29 and lattice methodJSOl, which indicate a gluon 
confinement mechanism equivalent to that in Landau gauge. Investigations in other 
gauges do not yet lead to a satisfactory understanding of gluon confinement even 



a The hermiticity assignment is valid in Landau gauge, although it is not the conventional on< 
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in the vacuum, as the assumptions made turned out to be too restrictive^^. Thus 
it is still unclear how the observations collected here translate to other gauges. 

A central issue is the fate of confinement and what "deconfinement" implies. 
Therefore adequate criteria are needed to test confinement. These are partly encoded 
in the pertinent 2-point functions. The releva nt p roperties will be only listed here. 
A brief review of them can be found elsewhere"^. 

Empirically, if the spectral function of a particle is not positive semi-definite, 
no Kallen-Lehmann re pres entation exists. It is then not part of the physical spec- 
trum and thus confined 3 - 3 -! A violation of positivity occurs, if the corresponding 
propagator D vanishes at zero momentum 

lim D(p 2 ) = 0. (2) 

p 2 ^0 

This is a sufficient, but not a necessary condition. Intuitively, for a massless particle 
condition (J2J yields a vanishing would-be on-shell propagator. Thus the correspond- 
ing particle does not propagate and is confined. 

The behavior in cq. © is also predicted by two confinement scenarios. The one 
of Kugo and Ojima° puts forward the idea that all colored objects form BRST 
quartets and are thus unphysical. One precondition for this is an unbroken global 
color charge. In the Landau gauge, this condition can be cast intoSS 

lim p 2 D G (p 2 ) -> 00, (3) 

where Dq is the propagator of the Faddeev-Popov ghost. Such an infrared diver- 
gence relates to long-ranged spatial correlations, stronger than those induced by a 
Coulomb force, since the divergence in momentum space is stronger than that of a 
massless particle. In this sense, ghos ts media te confinement. 

In the Gribov-Zwanziger scenaric0^EII entropy arguments suggest the domi- 
nance of field configurations close to or on the Gribov horizon in field configuration 
space. The Gribov horizon in Landau gauge is characterized by zeros of the inverse 
ghost propagator, predicting eq. ©. For an infrared constant ghost-gluon vertex, 
which is supported by lattice simulatio ns!^ and calculations based on 
eq. (01 follows for the gluon propag ator 23 . A further consequence of this scenario 
is that the truncation introduced below becomes exact in the infrared. 



3. Finite temperature and Dyson-Schwinger equations 

The basic quantity encoding the content of a theory is the infinite set of Green's 
functions. Their knowledge allows to calculate any thermodynamic quantity, as they 
completely determine the partition function. 

In practice, however, only a finite number of Green's functions can be determined 
within approximation schemes. These can be used to estimate thermodynamic quan- 
tities using e.g. the Luttinger- Ward/ Cornwall- Jackiw-Tomboulis formalism^! The 
thermodynamic potential can alternatively be determined by lattice calculations. 
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However, in most cases its microscopic origin remains elusive. Furthermore, as in- 
frared singularities are anticipated according to condition J3J, a complementary 
continuum formulation, not plagued by finite volume effects, is desirable. 

The Green's functions can be determin ed by the DSEs, which are the quantum 
analog of the classical equations of motio n * I Finite temperature can be intro- 
duced into the equations using the Matsubara formalism^! which compactifies the 
Euclidean time dimension to a length of 1/T. Thus, the zero component of the 
momentum becomes discrete with po = 2TrnT for bosons, where n is an integer in- 
cluding zero b . Furthermore, Lorentz invariance is no longe r m anifest^! Therefore, 
the propagators of the ghost and gluon can be described a 

Da(p) = _°6*£> (4) 

DM = PrM^fi 1 + PlM^^ 1 • (5) 

Here Pl and Pr are projectors longitudinal and transverse to the heat b ath, re spec- 
tively. Both are 4d-transverse, as required by Slavnov-Taylor identitie a ' \ The 
soft mode (po = 0) of Pt consists only of space-space (chromomagnetic) components 
while the one of Pl only contains time-time (chromoelectric) components. G, Z and 
H are scalar, dimensionless dressing functions, fully specifying the propagators. At 
zero temperature, Z and H coincide. 

The DSEs form an infinite set of coupled, non-linear integral equations for the 
Green's functions. By expanding the Green's functions in powers of g ordinary 
perturbation theory is recovered. If any progress is to be made beyond resummcd 
perturbation theory, it is necessary to solve these equations self-consistently. This 
is in general not possible, ft is therefore necessary to truncate the system to 
make it tractable. Truncations are in general a-priori not justifiable, ft is there- 
fore necessary to check their validity a -posteriori. The truncation scheme used 
for the results presented has been validated in various systematic 

wav JiW^ I ^ I 4 ^ 4 ^ 4 ^ Some of the assumptions have been tested using 

lattice calculation J'^ l ^ 1 1521 l n addition, the resul t s are in ve ry good agreement both 
qualitative and quantitative with available latticJ^ESEEl anc j calculationJ^EI. 
The truncated equations are graphically represented in figure ^ 

A serious ob stacle in any non-perturbative calculation for gauge theories is the 
Gribov problerrP^E^, as a local gauge condition is not able to fully specify the 
gauge. However, recently is has been argued^ that it is sufficient to stay within 
the first Gribov horizon to evade this problem. This is achieved by requiring the 
dressing functions G, Z, and H to be positive semi-definite. The influence of Gribov 
copies has been studied on the lattice and no contradiction to this assumption has 
been found^l 



b Note that this also applies to ghosts d espi te their Grassmannian nature, as they live on the 
periodic representation of the gauge groupS2J 
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Fig. 1. The truncated Dyson— Schwinger equations at finite temperature. The dotted lines denote 
ghosts, the dashed lines longitudinal gluons and the wiggly lines transverse gluons. Lines with 
a full dot represent self— consistent propagators and small dots indicate bare vertices. The open 
circled vertices are full and must be constructed in a given truncation scheme. A bare ghost— gluon 
vertex and slightly modified bare gluon vertices have been used here. Note that for soft modes the 
ghost-longitudinal and the 3-point coupling of three longitudinal and of one longitudinal and two 
transverse gluons vanish. 



A further interesting property is obtained in the limit of infinite temperature. 
This is equivalent to the static limit of the theory, as the exte nsio n of the time direc- 
tion shrinks to zero. Thus it can be described by a 3d theory^. The corresponding 
3d Lagrangian can be obtained as an effective field theory**!. In this course the orig- 
inal Aq component of the gluon field becomes an additional adjoint Higgs field, thus 
preserving the number of degrees of freedom, as a 3d gluon field has only one trans- 
verse polarization. The constants that appear in the 3d theory , like the Higgs mass, 
can be obtained by matching with the original theorjEHESES The most important 
quantity in this context is the (now dimensionful) 3d gauge coupling 173. Within 
the truncation scheme employed here™, it is given by g\ = g 2 T. The value of the 
running coupling g depends on the renormalization scale /i that can be selected at 
will. Usually fj, is chosen to be proportional to T, yielding at high temperatures 
a logarithmically decreasing coupling due to asymptotic freedom. However, to ob- 
tain a well-defined infinite-temperature limit wi th a fi nite 3d gauge coupling 173, it 
is here necessary to fix the value of g T insteaoPEl 

Thus the running coupling 

decreases polynomially with temperature. 



4. Results 

4.1. Asymptotics 

For asymptotically small and large (Euclidean) momenta, the truncated DSEs can 
be solved analytically. In the ultraviolet this is due to asymptotic freedom. Thus 
at momenta large compared to Aqcd an d T. th e same logarithmic running of the 
dressing functions occurs as in the vacuu m ™ Z and H coincide in this regime . 
Using an appropriate three-gluon vertex, resummed perturbation theory to one- 

c At sufficiently small temperatures this is already the case in the non-perturbative domain. 
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loop order is also recovered in the self-consistent solutions of the DSeJ^SI. 

This behavior persists at all temperatures. At sufficiently large temperatures, 
T Aqcd, the intermediate momentum regime Aqcd < P < T, becomes indis- 
tinguishable from the perturbation theory of a 3d theory of the soft modes alone. 
This is due to the decoupling of the hard modes, which have an effective mass of 
at least 2irT. This is a hard scale. Thus the hard modes are essentially perturba- 
tive and the Appelquist-Carrazone theorem applies. When the infinite-temperature 
limit is taken, this affects the full momentum range, and the effective 3d theory is 
obtained^! This behavior is explicitly seen in the DSE resultJZDHI 

In the infra red, the tr uncated system of DSEs can be solved analytically and 
self-consistentlySEHEIIlZI 

by power-law ansatze for the soft modes, 

G(p) = A g (p 2 )-« (6) 

Z(p) = A z (p 2 )- t (7) 
H(p) = A h {p 2 )- 1 . (8) 

The exponents are related^. In the vacuum —2k — t and t = I, and it can be 
shown that k > quite generally^! However, the exact value of k depends on the 
truncation^!, but strong arguments exist for k > 1/2. Therefore conditions J5J and 
® are satisfied and gluon s are confined in the vacuum. The best value to date, 
in agreement between DSlJSZEl anc j RG methodJ^, is k « 0.595, a value which 
is compatible with latest lattice resultJ^, suggesting k rj 0.53. Numerical results 
from DSE calculations at non-zero temperatures below the phase-transition find 
no change of the exponential. 

Interestingly enough, the terms stemming from the gauge-fixing part of the 
Lagrangian l|T]l. i.e. only diagrams on the r. h. s. in figure ^with a t least one ghost 
line, dominate in the infrared. This agrees with the Gribov-Zwanziger scenario and 
thus a Yang-Mills theory would be a topological field theory of Schwarz type in the 
infrared. Therefore no propagating modes exist and confinement is manifest. 

The situation in the high-temperature phase is different. The value of the expo- 
nents are again independent of the tempe ratu re and can most directly be calculated 
in the dimensionally reduced theorjE3"ZGID j n this case —2k = t + 1/2. Two solu- 
tions are found with k = 1/2 and k « 0.39, both satisfying conditions and 
The phase is therefore strongly interacting and confining, even in the limit of infinite 
temperature, as was anticipated^! ConsiderationJ2HE3 concerning smoothness in 
dimensionality and indications from the thermodynamic potential favor the value 
K w 0.39 as do lattice result J^. Again only the gauge-fixing part contributes in 
the transverse sector. However, in the longitudinal sector the situation is different. 
Here t ^ I due to interaction o f the so ft modes with the hard modes, leading to the 
generation of a screening masa23E2l A massive particle has an infrared cons tant 
propagator, requiring I = — 1. This result also agrees with lattice calculation^^!. 
Thus, longitudinal modes in contrast to transverse modes are not confined accord- 
ing to condition J2Jl. Therefore the high-temperature phase is qualitatively different. 
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Confinement in the longitudinal sector is discussed further in section T4. 31 

The hard modes, on the other hand, are restricted to p 2 = Po+p 2 > Po- Infrared 
is therefore measured with respect to p 2 . In this sense, the hard modes exhibit in 
both phases the behavior of massive particles, i.e. k% = th = lh = — 1 as expected. 
At sufficiently low temperatures non-perturbative effects still contribute, and the 
hard modes cannot be described accurately by perturbation theory alone. 

An interesting feature of the results 1)618(1 is a qualitative independence of the 
gauge group: For any semi-simple compact Lie-group, such as the physical SU(3), 
the results are qualitatively the same. This is most pronounced in the infinite 
temperature limit, where the results, after proper rescaling, are also quantitatively 
independent of the gauge group. This property hinges on the validity of the approxi- 
mations made, as here the gauge group only enters through its adjoint Casimir. This 
factor drops out (as does the gauge coupling) in the determination of the infrared 
exponents. It is as yet unclear whether this property persist beyond the present 
truncation. Therefore, in the following only the gauge group SU(3) will be consid- 
ered, except for the comparison with lattice data, using SU(2) where necessary. If, 
however, the truncation made here is indeed exact in the infrared, as is suggested 
by the Gribov-Zwanziger scenario, all such gauge groups have a qualitative unique 
infrared limit. This also implies that the mechanism of gluon confinement and the 
phase structure in Landau gauge is independent of the gauge group, which is an 
interesting observation to be investigated further. 



4.2. Full solutions 

The connection between the analytical infrared and ultraviolet regime has been 
obtained by lattice calculations as well as by self-consistently solving the DSEs 
numerically. There are various techniques for the latter case. In the low-temperature 
regime a torus discretization has b een use d ^ l ^'-'^ . In the high-temperatur e regime 
a continuum method was employecP^E^I. In the vacuum both have been usecPE"^l 

The results for the soft modes as a function of temperature are shown in figures 
121 II For the h igh-te mperature phase only the more likely solution is shown. The 
other is similai ^ , even quantitatively. The hard modes are essentially tree-level 
up to small corrections and are t hus not show n here. Results for these can be found 
elsewhere for temperatures below^l and above^EEl th e phase transition. 

In figure[21the ghost dressing function is displayed. The divergence in the infrared 
is clearly visible, satisfying Q. The dressing function does not exhibit a marked 
qualitative dependence on temperature. Thus, the ghost is able to mediate long- 
range forces at all temperatures, in agreement with available lattice calculations. 

The transverse gluon propagator is shown in figure |21 It is infrared suppressed 
at all temperatures, satisfying J3J). Therefore transversely polarized gluons are al- 
ways confined. Hence, the Yang-Mills theory remains non-trivial in the infinite 
temperature limit. The peak at mid-momenta roughly signals the transition from 
perturbative to non-perturbative behavior. The height and to some extent the posi- 
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Fig. 2. The ghost dressing function is shown in the vacuunB land at temperatures below the phase 
transitiorl^l in the left panel and above the gha.se transitiorESEU in the right panel. The lattice 
data in the left panel are at zero temperatur e!^ ! The three momentum is signified by P3 = \p\- 
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Fig. 3. The transverse gluon propagator is shown in the vacuu nil and at temperatures below the 
phase transitior!22J in the left panel and above the phase transitiorGSEI] in the right panel. The 
lattice data in the left panel are at zero temperature from a 20 3 X 64 latticed and in the right 
panel at infinite temperature continuum extrapolated^! and from large latticeJ^l 



tion of the peak is the region most sensitive to the truncation. The low-temperature 
propagators in both approaches do not fully reach into the infrared due to numeri- 
cal limitations. In addition, the lattice data in the vacuum do not bend over due to 
finite volume effects. Only recent calculations on an extremely asymmetric lattice 
of 16 3 x 256 find such a bending over®. Thus further investigations are of great 
interest. In the infinite-temperature limit, advantage can be taken of the lower di- 
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Fig. 4. The longitudinal gluon propagator is shown in the vacuu and at temperatures below 
the phase transitiorl^l in the left panel and above the phase transitiorESEU in the right panel. The 
latt ice data in the left panel are at zero temperaturJZland in the right panel at infinite temperature 



mensionality, and much larger lattices can be used. Such calculations show a clear 
bending over in the infrared^. 

The situation is different for the longitudinal gluon propagator as shown in 
figure^] The low-temperature phase is confined, and the longitudinal and transverse 
gluon propagators are similar. Here, also the DSE results are not suppressed in 
the infrared. This is likely due to finite volume effects in the numerical method, 
as close investigations sho w ^9 1 The weaker suppression in the infrared is more 
likely due to an increase of the coefficient Ah in (JHJ rather than due to a change 
of the exponent /. More important is the behavior in the high-temperature phase. 
The propagator goes to a finite value in the infrared, in agreement with lattice 
results. Thus the propagator of gluons polarized longitudinally w.r.t. the heat-bath 
is essentially dominated by a (renormalized) mass. 

The main difference between both phases is thus the longitudinal sector. This 
agrees with investigations of the temporal and spatial string tension"^. The former 
corresponds to the longitudinal sector and vanishes in the high-temperature phase. 
The spatial string tension corresponds to the transverse sector and even inc reas es. 
This also complies with a broken -/^-symmetry in the high-temperature phasJ^, as 
the corresponding order parameter is only observe d usin g time-like Polyakov lines. 
Furthermore, the lattic c 26 62 1 and the DSE resultsPSESI suggests that the infinite - 
temperature limit is essentially reached already at temperatures of a few times T c . 

At this point the results can be compared to the expectations. Originally, the 
generation of chromoclectric and chromomagnetic screening masses has been ex- 
pected. Long ago it was already showrl^that a magnetic screening-mass can be at 
best of order g 2 T and not of perturbative origin, if it exists at all. This problem lead 
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to the approach using which is based on resumming the hard mode contri- 

butions in self-energy diagrams, and then perform perturbation theory using these 
hard thermal loops. In the transverse infrared sector it is still plagued by severe 
problems, due to its perturbative nature, very similar to the vacuum case. In the 
high-temperature phase the HTL approach results in a transverse gluon propagator 
with a particle-like pole at p = 0, thus t — 0. 

The DSE results disagree with t hese fin dings and yield that a magnetic screening 
mass does not exist and is not likely^EHI. Instead over-screening persists as in the 
vacuum. This is in agreement with lattice calculations^! at the highest available 
temperatures. HTLs and the DSEs find qualitatively similar results only for the 
propagators of the soft longitudi nal mo de and the hard modes. However, the mass 
of the longitudinal mode is found23E5lto be of order g 2 T rather than of order gT, 
as expected in perturbative calculations. 



4.3. Positivity and confinement 

A sufficient criterion for confinement is the violation of positivity, which is already 
implied by the fulfillment of condition J5J. A more direct test is the determination 
of the corresponding Schwinger function, the Fourier transform of the propagators^. 
It is not positive semi-definite if positivity is violated. This has been found for the 
gluon propagator in th e vacu um^ and for the transverse gluo n prop agator in the 
high-temperature phasd^BEll j n agreement with lattice resultd^^™. 

As the longitudinal gluon propagator does not fulfill J5J in the high-temperature 
phase, it is of particular interes t to d etermine its Schwinger function. It is found 
to be not positive semi- definite™ 

This result depends, however, much more 
strongly on the truncation than in the transverse sector. Thus, further investiga- 
tions are mandatory to decide whether all gluons remain confined. The longitudinal 
gluon cannot be described by leading-order perturbation theory alone, and non- 
perturbative or higher-order pert urbativ e effects still contribute, as a comparison 
of lattice^ and DSE results shows^El. 

If the longitudinally polarized gluon is confined as well, this seems to occur 
through a different mechanism than for the transverse ly-p olarized gluon. The prop- 
erties of the Schwinger function indicate the possibilitjEBT f or a Gribov-Stingl-likJ"^ 
confinement mechanism, but the results are not conclusive. 



4.4. A note on the thermodynamic potential 

Lattice calculations find a thermodynami c poten tial which seems to logarithmically 
approach a Stefan-Boltzmann behavior I I I, but still differs significantly from 
it even at several times T c . Besides lattice calculations, several attempts have been 
made to obtain the thermodynamic potential. Perturbation theory yields a slowly 
converging expansion of t he potential, which differs from the ideal gas value even 
at very large temperatures^! An alternative is based on perturbation theory in the 
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hard and longitudinal modes and adding a few quantities from lattice simulations 
of the 3d Yang-Mills theory. This imp rove d g or weak-coupling expansion provides 
a good description of the lattice resultJ^. 

Using DSEs, the thermodynamic potential due to gluons can be obtained us- 
ing the Luttinger- Ward/ Cornwall- J ackiw-Tomboulis formalism^! It is, however, 
strongly limited by truncation artifacts. In the low-temperature phase, the thermo- 
dynamic potenti al is essentially constant, thus exhibiting the behavior expected of 
a confined systenP^. At sufficiently high temperatures, the potential is dominated 
by the hard mod es. These provide, up to small corrections, the potential of an ideal 
gas of gluonJ23, as the contribution due to the strong interactions is sub-leading. 
Hence despite the residual strong and confining interactions, a nearly trivial thermo- 
dynamic behavior is obtained. Near the ph ase t ransition, however, the contribution 
of the soft modes are likely to be import antP^I 

In principle, knowledge of the thermodynamic potential gives access to the criti- 
cal temperature, which is known from lattice calculation^ with high accuracy to be 
T c = 269 ± 1 MeV. Using DSEs it is not yet possible to extract the critical temper- 
ature. Nonet heless, observations which can be interpreted as super-heating^3 and 
super-cooling23ESI together with the qualitative change in the longitudinal sector 
suggest a first-order phase transition. This would be correct for SU(N) with N > 2, 
but not for other gauge groupJ^, e.g. SU(2). 



5. Discussion and conclusions 

Before discussing the implications of the results, it is important to assess their 
reliability. In the ultraviolet, resummed perturbation theory provides a reliable de- 
scription due to asymptotic freedom. Thus, the incorporation of the results of per- 
turbation theory is required for a final understanding. The infrared regime does 
not permit such a systematic approach. However, substantial evidence exists for the 
validity of the truncation scheme used in the DSE analysis "*^"^"* ™". These re- 
sults are reinforced by the comparison with available calculations using lattice and 
RG methods. The gap at intermediate momenta is well under control using lattice 
calculations, which agree qualitatively with DSE calculations. Thus the presence 
of gluon confinement in the vacuum and in the infinite -temperature limit is a ro- 
bust statement. This necessarily implies the non-triviality of the high-temperature 
phase of Yang-Mills theory, as expected. 

At finite temperatures, and especially in the domain of the phase transition, the 
results for the propagators become less reliable. Both, DSE and lattice methods, 
have problems to systematically treat this range. For DSEs, the problem is present 
mainly due to the necessity to handle an infinite number of Matsubara frequencies, 
while lattice calculations suffer from the required volumes to reach small tempera- 
tures. This also entails that the extension in the time direction becomes too small to 
reliably extract the longitudinal propagator. Nonetheless, the lattice provides excel- 
lent results on thermodynamic properties. On the other hand, some general results 
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on the infrared regime obtained from DSE calculations still hold. Thus qualitative 
results seem to be trustworthy. This is mainly possible due to the large effective 
'mass' of the hard modes, which is at the phase transition already of the order of 
1.7 GeV. 

Still, much work has to be done to obtain the quantitative properties, including 
the structure of the phases. In addition, it is as yet unclear which are the cor- 
rect colorless and thus gauge-invariant excitations governing the high-temperature 
phase. The persistence of confinement requires that this cannot be simple gluons, 
but perhaps light glueballs formed of soft gluons. These would likely not contribute 
significantly to thermodynamic quantities at sufficiently high temperatures. How 
this can be realized is subject to further investigations. 

Nonetheless, a coherent picture emerges. The main difference between the 
low-temperature and the high-temperature phase is not primarily one between 
a strongly interacting and confining system and one with only quasi-free parti- 
cles. Over-screening remains in at least part of the spectrum and thus confinement 
persists. In the transverse sector, the results are in accordance with the Gribov- 
Zwanziger and/or the Kugo-Ojima scenarios. The scenario that emerges is a chro- 
moelectric phase transition of first order. In the vicinity of the phase transition, the 
non-perturbative effects are likely relevant to thermodynamic propertied^, under- 
lining the importance for experiments. 

It is still unknown which eff ects can be induced by quarks. In the infinite- 
temperature limit, arguments exisl^that they will not affect the propagators due to 
their fermionic nature. Fermions require anti-periodic boundary conditional^, and 
therefore generate effective 'masses' of (2n + l)irT. Thus all quarks become infinitely 
massive in the infinite temperature limit, irrespective of their intrinsic mass and 
thus decouple due to asymptotic freedom. At the phase transition, however, lattice 
calculations indicatd-^ that the quarks possibly are able to change the nature of 
the phase transition to a second order one or to a cross-over. 

In addition, the observation of a drastic change in the inter-quark potential^ 
and its connection to the residual confinement of gluons still needs to be under- 
stood. It cannot as yet be firmly concluded that the inter-quark potential still rises 
slowly at very large distances even at high temperatures. On the other hand, quark 
confinement is still not yet understood at zero temperature^, thus posing a more 
significant problem. 

In this review the main results from studies on the finite-temperature behavior 
of Yang-Mills theories in Landau gauge have been summarized. The investigation 
is based primarily on the infrared properties of the propagators of the elementary 
degrees of freedom, the gluons and the Faddeev-Popov ghosts. 

The results provide evidence for the existence of a phase transition in the pure 
gauge sector. In both phases strong interactions and non-perturbative features 
like confinement are found. Direct comparison to lattice and RG results in the 
vacuum and at high temperatures support this. However, many unsettled ques- 
tions remain, such as the nature of the observable degrees of freedom in the high- 
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temperature phase and the influence of quarks. Thus further investigations are de- 
sirable and i n prog ress. These are accompanied by corresponding investigations at 
finite densit y '. This embeds these findings in the ongoing effort not only to 
understand the vacuum structure of QCD but also to map out the phase diagram 
and to finally understand the results obtained in experiment from first principles. 
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